Dephasing of quantum bits by a quasi-static mesoscopic environment* 
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We examine coherent processes in a two-state quantum system that is strongly coupled to a 
mesoscopic spin bath and weakly coupled to other environmental degrees of freedom. Our analysis 
is specifically aimed at understanding the quantum dynamics of solid-state quantum bits such as 
electron spins in semiconductor structures and superconducting islands. The role of mesoscopic 
degrees of freedom with long correlation times (local degrees of freedom such as nuclear spins and 
charge traps) in qubit-related dephasing is discussed in terms of a quasi-static bath. A mathemat- 
ical framework simultaneously describing coupling to the slow mesoscopic bath and a Markovian 
environment is developed and the dephasing and decoherence properties of the total system are 
investigated. The model is applied to several specific examples with direct relevance to current ex- 
periments. Comparisons to experiments suggests that such quasi-static degrees of freedom play an 
important role in current qubit implementations. Several methods of mitigating the bath-induced 
error are considered. 

* Dedicated to Anton Zeilinger, whose work has inspired exploration of quantum phenomenon in 
many avenues of physics and beyond. 
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I. INTRODUCTION 

Solid-state quantum information research attempts the 
difficult task of separating a few local degrees of freedom 
from a strongly coupled environment. This necessitates 
a clever choice of logical basis to minimize the dominant 
couplings and environmental preparation through cool- 
ing. Electron spin in quantum dots has been suggested 
as a quantum bilià, where orbital degrees of freedom 
do not carry quantum information. As a result, the in- 
dividual quantum bits are relatively well isolated from 
non-local degrees of freedom, interacting strongly with 
only a local spin environmen1jà£iàíi£i2i2ài2i and weakly to 
phonons through spin-orbit couplingàü*i2iüi4. Joseph- 
son junctions and cooper-pair boxes in superconduct- 
ing systems can be engineered for decreased coupling 
to local degrees of freedomià, and several different ap- 
proaches to superconductor-based quantum bit devices 
have shown remarkable promise as quantum computa- 
tion devicesiàiiíààiLià. However, superconductor designs 
may strongly interact with local charge traps and mag- 
nètic impuritiesi2i2Sí2i ; leading to errors from local noise. 

In this paper we examine the problem of coupling to 
the local environment for a solid-state system. By fo- 
cusing on quantum information-related tasks, the detri- 
mental effect of local degrees of freedom on specific op- 
erations can be assessed. In our model, so-called non- 
local degrees of freedomSà (such as phonons and pho- 
tons) are assumed to be weakly coupled and considered 
in the normal Markovian limits. However, local coupled 
systems corresponding to spins, charge-traps, and other 
finite-level systems, can be interpreted in terms of a fi- 
nite set of nearby spins that may be strongly coupled to 
the qubit. In some cases, such as an electron spin in a 
quantum dot, this heuristic picture corresponds exactly 
to the actual hyperfine coupling between electron spin 
and lattice nuclear spins. In those systems, the local en- 
vironment consists of 10 4 -10 6 lattice nuclear spins, and 
the coupling can be strongíüàiSà. For superconductor- 



based designs, the quantum bit couples both to actual 
spins and to charge degrees of freedom such as 1/ /-type 
fluctuat or s2», which can be modeled as two-level systems. 
Thus, analysis of the most general coupling of a qubit to 
a mesoscopic collection of spins yields an understanding 
both of limitations in current experiments and of meth- 
ods for improving the operation and design of solid-state 
quantum computation devices. We remark that several 
works have analyzed this type of situation in detail for 
specific implementations, in superconducting qubits22i2I 
and in quantum dot qubits2*i2iSàí2Si2i. We also note that 
recent work focusing on quantum dot-based systems par- 
allels several significant elements of the current paper—; 
the results were arrived at indcpcndcntly. The present 
work focuses on the generality of the model and its im- 
pact in the context of quantum information processing. 

We begin by considering a spin- 1/2 system (the 
"qubit") coupled to a finite number of spin- 1/2 degrees 
of freedom (the "mesoscopic bath") to both weak and 
strong driving fields. Inclusion of weaker, non-local cou- 
plings (the "Markovian environement" ) through a Born- 
Markov approximation reveals a natural hierarchy. A 
separation of time scales indicates that the local environ- 
ment is quasi-static with respect to experimental time 
scales with long-lived correlation functions. However, 
over many experimental runs, the local environment fluc- 
tuates sufficiently to defy efficient characterization and 
correction by means of direct measurement. The genèric 
nature of the coupling allows bath characterization using 
only a few parameters. Then, detailed analysis of spe- 
cific operations (phase evolution, Rabi oscillations, spin- 
echo) and comparison to experimental results is possible. 
Whilc such a spin-bath model has been considered previ- 
ous1jí2^Sí*22íS2éSÍ, the significance of the mesoscopic and 
quasi-static limits has not been emphasizcd. 

We find that a quasi-static, mesoscopic environments 
results in fast dephasing, induces power-law decay and 
non-trivial phase shifts for driven (Rabi) oscillations of 
the qubit, and may be corrected for using echo tech- 



niques. The generality of the model suggests it may 
be appropriate for many systems with strong coupling 
to a stable local environment. In addition, comparison 
to experiments suggest this model can help explain the 
reduced contrast and long coherence times of current 
solid-state qubit systems, indicating that the local en- 
vironments of current solid-state qubits both have long 
correlation times and play a crucial role in the dynamics 
of the qubit. Several methods of mitigating the effects of 
the bath are considered, all of which take advantage of 
the long correlation time to reduce errors. Some of these 
tcchniques are extensions of established passive error cor- 
rectior. 31 and quantum bang-bang ideas 3 ^, while another, 
the preparation of coherences in the environment to re- 
duce dephasing (environmental squeezing) is introduccd 
in this context for the first time. 



II. THE PHYSICAL SYSTEM 

Consider a relatively well-isolated two-level system (a 
qubit) in a solid-state environment. This could be a 
true spin-1/2 system such an electron spin or an anhar- 
monic system such as a hybrid flux-charge superconduc- 
tor qubits^ 3 ». The qubit is coupled to a collection of spin- 
1/2 systems comprising the mesoscopic bath. In addition 
to these potentially strongly coupled degrees of freedom, 
other, weaker couplings to bosonic fields (phonons, pho- 
tons, and large collections of spins that can be mapped 
to the spin-boson modeli) are included as a Markovian 
environment, and act as a thermal reservoir. Qualita- 
tively, this leads to a heirarchy of couplings. The qubit 
plus mesoscopic bath are treated quantum mechanically, 
while the additional coupling to the larger environment 
(the thermal reservoir) is included via a Born-Markov ap- 
proximation. Thus, the reservoir plays the role of ther- 
malizing the mesoscopic bath over long time scales and 
providing additional decay and dephasing of the system. 
By assuming slow infernal dynamics and thermalization 
for the mesoscopic bath, a quasi-static regime is investi- 
gated, in contrast with the usual Markovian approxima- 
tions for baths. 

Starting from first principies, the hamiltonian is intro- 
duced. A transformation to the rotating frame allows for 
adiabatic elimination which simplifics the interaction, af- 
ter which the quasi-static assumption is considered. An 
explicit bath model is chosen, and the role of the thermal 
reservoir included through a Linblad-form Louivillian. 

A. Hamiltonian 

The qubit, 2-lcvel system, is described by the spin-1/2 
operators S x ,S y , and S z . An external biasing field pro- 
duces a static energy difference between the two levels of 
angular frequency lo. This serves as a convenient defi- 
nition of z-axis. A driving field, with Rabi frequency Í2 
and oscillating with angular frequency v, is applied along 



a transverse axis in the x — y plane. In most systems, u> 
is fixed, and modulation of Q turns on and off rotations 
around the transverse axis. The relative phase of the 
driving field to the biasing field controls the axis about 
which rotations occur. We choose the phase such that 
rotations occur around the a;-axis. The hamiltonian is 

H sys = uS z + 2f2cos(z/í)«S X) (1) 

in units where h — 1. This hamiltonian is commonly 
encountered in quantum òptics, NMR, ESR, and solid- 
state quantum information devices. For solid-state qubit 
systems, u> and v are typically ~ 1 — 100 ns _1 and íl ~ 
0.01 - 1 ns _1 . 

Adding the most general environmental coupling pos- 
sible, the total hamiltonian is 

Htot — H sys + XS ■ E + He (2) 

where A is a coupling constant to the environment and E 
is a vector of environmental operators. The tilde terms 
indicate that this coupling includes both the mesoscopic 
spin bath and the Markovian environment. All of the 
internal environmental dynamics, including thermaliza- 
tion, is encompassed in He- 

Separating the environment E into the mesoscopic spin 
bath and Markovian environment, the interaction and 
environment terms in the hamiltonian can be rewritten: 

XS-É + He = XS ■ À+ X C S ■ C + H A + H c + È AC . (3) 

The A terms denote the quasi-static bath degrees of free- 
dom, while the C terms correspond to a larger environ- 
ment that can be treated in a Born-Markov approxima- 
tion. Qualitatively, the complete picture of bath and en- 
vironment drawn in this paper is hierarchical in nature 
(Fig.[TJ, with the system (qubit) coupled to a mesoscopic 
spin bath, both of which are coupled to larger environ- 
ments. 

This paper focuses on the regime of strong system- 
bath coupling and weak system-environment and 
environment-bath coupling. In this limit, we can treat 
the Markovian environment 's effects entirely through 
a Lindblad form Louivillian acting on the combincd 
system-bath. We write the effective hamiltonian as 

H eif =È ays + \g-ï+H A , (4) 

while superoperator describing evolution of the system- 
bath density matrix p is given by the differential equation 

p = i\p, Hes] + CcnvP ■ (5) 

The Louivillian, £ onv , is described below. 

Regarding other limits, the weak system-bath limit can 
be treated pertubatively when the bath-environment cou- 
pling and/or system-environment coupling is stronger. 
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FIG. 1: Heirarcichal system-bath-environment coupling. 
External control of the system (qubit) and measurement are 
independent of the environment The mesoscopic spin bath is 
strongly coupled to the system, while the environment pro- 
vides weak dephasing of the system and spin bath degrees of 
freedom. Infernal dynamics of the spin bath can be ergodic, 
but some experimental control of the bath may be achieved 
in limited situations (e.g., through NMR) as indicated. 



When the system-environment coupling can be experi- 
mentally controlled, such as during qubit manipulation, 
the system degree of freedom could be used to transfer 
cntropy from the bath to the environment in a controlled 
manner, yielding non-thermal initial states of the meso- 
scopic spin bath. The usefulncss of this final point will 
be discussed later. 

In many cases, the bias field u> is large (|u;| ^> |fl|, |A|). 
This is similar to the quantum optical case; accordingly, 
the Standard recipe of transformation into a rotating 
frame and a rotating wave approximation is also vàlid 
here, and off- axis (flip-flop) terms are removed, resulting 
in a Jaynes-Cummings-like hamiltonian. Roughly speak- 
ing, energy cannot be conserved when the system spin 
flips through changes of the bath state. However, the 
additional coupling to the bath leads to terms not found 
in the Jaynes-Cummings hamiltonian. 

More explicitly, applying the unitary transformation 
U = exp(— iu)S z t), the hamiltonian is 



XS Z À Z + \/4[e iut S + À_ 
Ha- 



'SLA, 



(6) 



When the detuning of the driving field is small, 5 = (v — 
lu) -C u>, the highly oscillatory terms in this rotating 
hamiltonian can be neglected. To do so, the propagator 

Ü(t) = Texp(— i J Q H(t')dt') is formally expanded over 
one period of the applied unitary transformation, r = 



4it/uj, and a Magnus expansioni gives 
Ü(t) 



cxp[-ir(^2 H k )] 

fc=0 



(7) 



Sl/2[S + +S-] + (5 + \A,)S s +H B (8) 

-i/2r f {\É{t 2 ),H{ti)] dt 2 dt x . (9) 
Jo Jo 



Ho = I/t H(t) dt 



Hi 



The first term Hq will be used hereafter. Neglecting H\ 
and higher terms of the expansion if formally equivalent 
to the rotating wave approximation used to derive the 
Jaynes-Cummings hamiltonian. 



B. Quasi-static bath assumptions 

In Standard analyses, usually the bath dynamics are 
fast relative to the system dynamics, and a Born-Markov 
approximation is appropriate. However, when the infer- 
nal bath dynamics are much slower than the system dy- 
namics, the infernal bath corrclation functions are long 
lived, and we can consider a quasi-static bath. In this pic- 
ture, one collective bath operator (here A z ) has, in the 
Heisenberg picture, a long-lived correlator. In particular, 



(À z (t + T)À z (t)) 



< 1 



(10) 



for t > 1/Í2. The quasi-static limit is approached when 
this is satisfied for t > í oxpt (the time for a single exper- 
imental run) but not for r ~ í t ot (the time to generate 
cnough runs to characterize A z ). 

We consider this in detail. Infernal dynamics of the 
spin bath and coupling to the Markovian environment 
will result in a decay of correlations. A general descrip- 
tion of the correlator is 



(À z (t)À z (t')) = / doSM^C-*') 



(11) 



where the spectral function S(üj) has some high frequeney 
cutoff, r, set by the infernal dynamics of the bath^&. 
When T 3> íl, the Markovian limit may be reached, and 
the quasi-static theory developed here is no longer nec- 
essary. 

When r — > 0, the thermalization time goes to infinity. 
Then, while the state of the bath is initially unknown, a 
series of measurements can be used to estimate the de- 
tuning of the system, and thus the value of XA Z (see, e.g., 
Refs. 25 37). If the manipulation process has a measure- 
ment time for a single projection of r m > 1/^, after n 
such measurements, the error of (A z ) will be 



A«A)) 



A 



(12) 



This procedure is exactly that of a classical phase esti- 
mat ion algorithm. It takes a time oc n to generate a \fn 
improvement in the knowledge of the current value of A z . 

In between lies the quasi-static limit, where a single 
experiment (i.e., preparation of a given qubit state, evo- 
lution under Hq for a specified time, followed by pro- 
jective measurement) can be performed well within the 
correlation time, but not enough experiments can be per- 
formed to well estimate A z . Qualitatively speaking, this 
is equivalent to taking a thermal average for the ensemble 
of measurements, and mathematically the same as the re- 
sult for the simultaneous measurement of a large number 
of exactly similar systems. The uncertain knowledge of 
A z leads to dephasing when attempting specific quantum 
information operations. 

As a typical case, consider an infernal bath hamil- 
tonian Ha — fi J\ ^B.k-^z ano - a coupling with 

A = a.kB k . This could describe nuclear spins for 
quantum-dot systems, or charge traps. As [íf^i, A+] = 
J~] k ctk^>B,kB^_, the rotating-wave approximation is vàlid 
if u>B,k <C w f° r a h k. Thcn, the application of a large 
bias field u>, e.g., due to a static magnètic field, makes 
the interaction essentially static and classical, whereby 

the system experiences a fluctuating (classical) field A; 
fluctuations are due to corrections to the rotating wave 
approximation and coupling to the Markovian environ- 
ment. Non-classical correlations of the bath, produced by 
off-resonant interactions»^ or preparation techniqucs 
(e.g^à) return the problem to the fully quantum domain. 

For example, for electron spin in quantum dots, where 
the bath is lattice nuclear spins, dipolar diffusion pro- 
cesses lead to flip-flop interactions wherein a spin in- 
side the dot is exchanged with that of a spin outside 
the dot, producing ergodic dynamics within the bath on 
ms timescales^Mü. In contrast, spin-lattice relaxation 
(coupling to the larger environment) can be minutes to 
hours in time. Over time, the dipole diffusion destroys 
the correlations of the A z parameter. Detailed calcula- 
tions, using the method of moments, suggest that dif- 
fusion in dots is within an order of magnitude of the 
bulk diffusion rate^i, and this diffusion decorrelates A z 
no faster than the time scale set by the NMR linewidth 
(10 ms -1 for GaAs^à). Thus the primary effect of spin- 
lattice relaxation (coupling to a Markovian environment) 
is to produce an initial thermal state of nuclear spins, 
while dipole-dipole interactions lead to decorrelation of 
A z (infernal bath dynamics). 

For concreteness, we choose for our mesoscopic spin 
bath of N spin- 1/2 systems to have the forní 



(13) 



where I z is the spin component corresponding to the kth 
bath spin, and an internal Hamiltonian giving the quasi- 
static limit. As many bath variables scale as \Na 2 , we 
normalize the au such that a 2 — Í/N, thus keeping A as 



the sole bath strength parameter. The bars above quan- 
tities denote statistical averages over static variables. As 
we have assumed all bath spins are spin- 1/2, we can 
choose the z axis without loss of generality. 



C. Markovian environment effects 

The principle effects of the Markovian environment are 
two-fold: it leads to relaxation of the qubit system, and 
it leads to long-time thermalization of the mesoscopic 
spin bath. We have established that the internal dynam- 
ics of the quasi-static bath, combined with coupling to 
the Markovian environment, will lead to decorrelation 
of A z and dynamic thermalization. For a weakly cou- 
pled Markovian environment, with equilibriation times 
for the bath that are much longer than experimental 
times, the main effect of the Markovian environment is 
to produce an initial mesoscopic spin bath state of the 
form pa oc e~ HA / ksT . 

In general, the additional Markovian environment de- 
grees of freedom will also couple to the system, leading 
to relaxation and decoherence that we can model with a 
Lindblad form of the Louivillian of the combined system- 
spin bath density matrix. Formally, tracing over the en- 
vironment yields a superoperator (see Eqn. [5J 



2 r y s [2S^pS+ — pS+S- — S+S-p] 



2I_pI + - pl + l_ - I+I_p 



+2I+pI- - pl-l+ - I-I+p 



(14) 



corresponding to the Markovian environment. We pre- 
sume that the spin-bath is in the high temperature 
limit, while the qubit system only has spontaneous de- 
cay. This is consistent with the rotating wave approxima- 
tion. The resulting equations, with pure radiative decay 
for the system, are cumbersome to work with, and we 
will sometimes use the input-output operator formalism, 
which is formally equivalent by the quantum fluctuation- 
dissipation theorem. 



III. FREE EVOLUTION 

In this section, free evolution for a variety of 
experimentally-relevant situations is considered. Free 
evolution corresponds to the case of = 0, i.e., the ab- 
sence of driving field. In quantum information, free evo- 
lution is equivalent to phase rotations (Z rotations), and 
the behavior of free evolution is an indicator of how good 
a quantum memory the chosen qubit provides. The decay 
of coherence during free evolution due to the mesoscopic 
spin bath may be severe; however, we consider spin-echo 
methods as a means of mitigating the bath's effects. 

Equation[H]is the bàsic Hamiltonian used hereafter; A z 
is given by Eqn. For free evolution, Í2 = 0, and the 



Hamiltonian is 



H 



FID 



(XÀ Z + S)S Z 



(15) 



For time-dependent À z such that [À z (t),À z (t')] = (ap- 
propriate for the specific baths consider in section II), the 
exact propagator from t = to time t is 



U FID (t) =exp[-i(St + X / À z (t')dt')S. 



(16) 



In the quasi static limit, wc may usually replace 
L À z (t')dt' with À z t, simplifying evaluation of the prop- 
agator. 



A. Free induction decay 

To illustrate free evolution, we start by studying free 
induction decay (FID) , considered for this case by2& and 
for a dynamical spin bath by^2>. FID corresponds to the 
decay experienced by a spin coherence (e.g., a density 
matrix element | t)(l I) during evolution in a Ramsey- 
type experiment. The system is prepared in the state 
|+) = 1/V2QT) + 14», with (S+(Q)) = 1/2, for example, 
by rotations induced by control of íl (considered below). 
The driving field is turned off (SI = 0) and the system 
evolves according to the S z coupling only. Thcn, 



4> 



FID 



(s+(t)) 



ex.p[-i(XÀ z + S)t] 



bath 



(17) 



For an uncorrelated bath ((Ï Z Ï 3 Z ) = {ï z ){ï j z ) for k ^ j) 
we can write this in terms of individual bath couplings 44 



FID 



i 61 



fi [cos(Àa fe í/2) - i{2Ï k z ) sin(Aa fe t/2)]. 



The minimum of the projection of the actually final state 
of the desired final state, Ffid = l 'E'í 1 / z> | 2 is bounded 
from below by (Mm k [\(2Í*)\]) N . Thus, initial bath po- 
larization for finite N limits the maximum decay. 

The initial decay of this coherence is quadratic, giving 



'FID 



(Aí) 2 



Y,^ k {l-{2Ï k z )% (19) 



as shown in the inset of Figure For a large number of 
spins, the intermediate time decay converges to a Gaus- 
sian, exp[— (7^/£ií) 2 /2], with a rate 



1pid = X /5> 2 fc (l-(2J^). 



(20) 



When the inhomogeneity in the ctf. coefficients is low, 
the system experiences mesoscopic revival on a time scale 




2 3 
t [/N/À] 



FIG. 2: FID evolution (F FID ) in time seales VN/X for 
increasing inhomogeneity, a a = {0.1, 0.3, 0.6, lj/V^V, = 
30, P = 0. The inset shows increasing decay in the 
short time limit due to increasing inhomogeneity, a a — 
{5, 10, 20} /VW,N = 30,P = 0. 



given by the single spin coupling strength, Aa ~ X/\/N . 
This type of revival is shown in Figure [3 where a Gaus- 
sian distribution of of increasing widths are used, and 
a bath polarization P = (2I Z ) = is assumed. These 
mesoscopic effects for finite spin systems have been in- 
vestigated in great detail elsewheraSà, and we refer read- 
ers to those works for a detailed description. In most 
physical settings, inhomogeneity is large and revival will 
not be observed. However, for systems where N is small, 
fluctuations can become substantial even with a large in- 
homogeneity in the ccfc's, and such fluctuations have been 
observed experimentally 4 ^. 

To put the results in an experimental context, we con- 
sider a Ramsey fringe experiment. A system is prepared 
in the |1) state. A 7r/2 pulse is applied along x (denoted 
iï^, 2 ), and then at a later time í, an opposite, —n/2 
pulse is applied. Afterwards, a projective measurement 
on z measures the rotation of the S y eigenstate due to en- 
vironmental effects and detuning. In terms of actual ex- 
perimental implementations, preparation of a spin-down 
eigenstate and measurement along the S z axis is enough, 
if we use two 7r/2-pulses with a delay t. Formally, in the 
limit of perfect pulses the propagator is 



U FID (t) = R x _„ /2 eM~^A z S z t]R 



71-/2- 



(21) 



In the limit of perfect pulses, the probability of measuring 
the final state in the | J.) state is given by Ffid as a 
function of the delay time, t. If the system is coupled to a 
spin bath, it will show Gaussian-type decay in the initial 
time period. This decay is limited in extent for baths 
of finite polarization (non-infinite temperature). For a 
small bath size, mesoscopic fluctations and revival after 
the initial decay may be evident. 

For a series of sequential measurements, the finite cor- 
relation time of the bath will be evident. For example, for 
fixed interaction time t, a series of measurements {Mk} 
giving 1 for | l) and for | f ) at times í& will manifest 



correlations due to the correlations of À z at the differ- 
ent times, tk- Assuming the quasi-static limit, such that 
M k = cos 2 (AÍ z (í fe )r/2), we find 

(M,M k )-(M) 2 = l e xp[-(Xr) 2 í duS{üo){e^-^ -l)] 

(22) 

where the cutoff T of S{ui) is such that IV <C 1- Thus, 
the correlation function for a series of measurements 
with rr <C 1 in principio allow reconstruction of the 
spectral function describing A z . We note that "/fi d = 

X^J J S{u))du). 



results of a spin echo experiment allow an alternative 
method of extracting the correlation function associated 
with the mesoscopic spin bath. For a spectral function 
describing the bath with cutoff T <C 1/í, the echo will 
decay according to exp[— 27Í — A 2 r 2 í 4 /4]. 

For free evolution (Í7 — 0), used for phase gates 
and quantum memory, we find that a mesoscopic spin 
bath leads to dephasing in the quasi-static limit, with a 
time scale given by 1/A. Sequential measurements in a 
Ramsey-type experiment allow for reconstruction of the 
of the spectral function associated with the bath. Fur- 
thermorc, spin-echo can be used to reduce the dephasing 
dramatically. 



B. Spin-echo decay 

As the bath correlation time is long, spin-echo tech- 
niques can exactly cancel this type of dephasing. For 
spin-echo, a 7r-pulse is applied in between the two 7r/2- 
pulses; the time delays before and after the 7r-pulse are 
t\ and t 2 , respectively. The total evolution is thus 

UsEfah) = K /2 exp[-i\À z S z t 2 ]K 

expHAÍAíi^/2 (23) 

exp[-iAÍ 2 S z ii]iÇ /2 (24) 
= R*_ n/2 expl-iXÀAih - t 2 )]Rl /a (25) 

which is the original FID propagator (Eqn. I21f) with 
t — > (íi — t 2 )/2. The previous results hold but now as a 
function of the time difference. In particular, no decay 
occurs when íi = t 2 \ The probability of measuring | j) is 

F SE (ti,t 2 ) = (cos 2 [AÍ,(í 1 -í 2 )/2]) 

= Í{l + exp[- 7 2 /í5 (í 1 -í 2 ) 2 /2]}.(26) 

Spin-echo fidelity (FsE{t, í)) is limited both by the im- 
perfections in the rotations, R^/ 2 and R%, Markovian 

processes that directly lead to decay of S z , and decay 
of the correlation function of A z , i.e., in the interaction 
picture, A z (t) ^ A z (t') for t much later than i'. The 
first type of error, examined in detail below, does not 
depend on the total time t\ + t 2 , while the latter errors 
do. Neglecting imperfections in rotations, the effect of a 
correlated bath on spin-echo is a textbook problem, and 
we merely cite the result here: 



F SE (t,t) = i 



,-2 7 t 



(27) 



xexp[-A 2 / S{uj)úxí A {utl2)/(uj/2fduj\ 



where we assume À z is a Gaussian variable described 
by the spectral function S(u>), as per Eauation llll The 



IV. DRIVEN-EVOLUTION 

While the free evolution of the system is illustrative of 
quantum memory, a critical operation for quantum infor- 
mation processing is rotations about an axis perpendicu- 
lar to the applied bias. This is achieved by a pulse in the 
transverse coupling field strength, O, leading to driven 
oscillations between ||) and ||). These Rabi oscillations 
form the basis for the x-axis rotations used, for exam- 
ple, in Ramsey-fringe experiments, and more generally 
for producing a universal set of single qubit gates. 

In practice, the perfect rotations (Rabi pulses) needed 
to produce perfect rotations are unavailable, due to the 
finite power availablc for such pulses. We now consider in 
detail rotations obtained by driving the system (non-zero 
fí in Eqn.0) in the presence of the mesoscopic spin bath 
and Markovian environment. In the Heisenberg picture, 
the single-spin (system) operator equations of motion are 
given by 



dt 



= i[H ,S] 



(28) 



For time independent À z (quasi-static limit) we can solve 

these equations exactly, noting that d/ dt S — Con x S, 
where 

(29) 
(30) 

This evolution corresponds to the optical Bloch equa- 
tions, giving driven rotations of the qubit Bloch vector 
on the Bloch sphere about an axis n at frequency Cu. The 
dependence these quantities upon the state of the bath 
leads to shot-to-shot fluctuations, which broadens ensem- 
ble or sequential measurements. 

Alternatively, we may write the propagator for Hq as 




Ü(t) = cos(ó)í/ 2 )ï - 2i sin(ó)í/ 2 )^ ' S 



(31) 



where we have taken the quasi-static limit, such that 



Returning to Eauations !29l·l3*Ül we calculate the rota- 
tion of an initial z eigenstate with (S z (t — 0)) = 1/2, 
i.e., determine the Rabi oscillations expected for the ini- 
tial state 1 1)- The field O leads to population and coher- 
ence oscillations. Considering the z and y components 
together, we define 



/ 

C(í) 



(XÀ Z + 5) 2 

(XÀ Z + s) 2 + n 2 



Q 2 exp(-iyj (XAz + S) 2 + íl 2 t) 

(\À Z + s) 2 + n 2 ' 



(32) 
(33) 



Then, (2S z (t)) = Ro[(/+C(í))] and (2S y (t)) = Im[(C(í)>]. 
The value of (/) is a measure of the maximum con- 
trast of population oscillations, and also gives the steady- 
state population difference and lincshape; this is approx- 
imately maximal for S = —X(A Z ). (C(í)) gives the oscil- 
latory part of coherence and population. 
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FIG. 3: Lineshapes for a bath with N = 10 and increasing 
inhomogeneity [a a — 0.1(green),0.5(red),l(blue)]. The addi- 
tional line is N = 20, o a = 0.5 (black). 



p(A). The integrals involved are solved in Appendix IaI 
We note that the result is 



A. Steady-state behavior 

For long times, the steady state behavior will be given 
by (/) , which measures the population difference between 
t and!. This gives the measureable response of the sys- 
tem to CW excitations of the system. A simple estimate 
of (/) is provided by writing A z as (A z ) + 5A Z , and as- 
suming ((a + À z )- 2 ) = ((a + (À z )) 2 + (AÀl))~ l ). Then, 



1 



l + n 2 /[(S + X(A z )) 2 + X 2 (AAl)} 



(34) 



This has a màxima at the mean bath detuning, X(A Z ), 
and as such this point would be the observed zero in 
detuning. Differences from this observed zero, denoted 
S, gives an absorption spectrum of 



</> = 1 - ~ 



Ü 2 



S 2 + X 2 (AA 2 Z ) + n 2 



(35) 



which bchaves as a Lorenztian with a linewidth given by 
the combined power broadening and bath broadening for 
Ü « A(AÍ2)i/2. 

The approximations used for this simple estimate 
break down for large bath strength. To find this behav- 
ior in the strong bath limit, whcrc higher order moments 
are important, we need not assume that the bath den- 
sity matrix is uncorrelated, as before, but rather that 
the bath density matrix is diagonal in an eigenbasis of 
the bath operator, A z à^ This would be appropriate for 
nuclei interacting with an electron spin in a quantum dot 
at finite magnètic field, or for charge traps capacitively 
coupled to the system. 

We can evaluate (/) by tracing over the bath den- 
sity matrix, replacing p bath = Y,\a P*a z | XA z )(XA z \ with 



p(-S) y / 2iríl/u 



(36) 



for the lineshape. This has assumed that u^ 1 is a fre- 
queney scale over which the bath density of states is 
relatively flat, and as such can be Taylor expanded to 
second order. This time scale is 



40 



Qp"(-S) 



(37) 



In other words, driven oscillations can measure the bath 
density of states as a function of detuning, convolved 
with a kernel of width u~ l . We plot several lineshapes 
inFig.OJ 

Essentially, we find that the steady-state behavior of 
the system indicates a maximum of population in the 
state opposite the initial state at a detuning 5 = —X(A Z ). 
In addition, in the weak field limit, the line shape pro- 
vides a sensitive measure of diagonal components of the 
bath's density matrix, so long as a calibration of Çl inde- 
pendent of the bath exists. 



B. Time-dependent behavior 

We now consider the actual, observable population and 
coherence oscillations induced by a driving (Rabi) field. 
Looking at the explicit time dependence, we consider the 
envelope of oscillations given by |£(í)|. The initial decay 
is quadratic, (1 — (/))[! — (7n i oí) 2 /2], and when renor- 



malized to give maximal contrast, it has a width 



7a o 




1 - (/> 



(38) 



(39) 



\ 



(\a z + s) 2 + n 2 



' (\À Z + s) 2 + n 2 , 



This suggests that measurements of such short time de- 
cay will measure second order moments of the bath op- 
erators. 

While the short time behavior is consistent with Gaus- 
sian decay, as would be expected from a non-Markovian, 
non-dissipative bath, the long time behavior is a power 
law. 

The integrals involved are solved in Appendix es- 
sentially we solve for (((t)), and find 



(C(*))«p(-í) 



,-ínt 



(0 



(40) 



where the time scale for phase shift, and phase shift angle, 
are given by 



r(t) = y/t 2 + u 2 

6{t) = 1/2 tan" 1 [t/u] , 



(41) 
(42) 



and u is the same as Eqn. 13*71 

We note that this result is much more robust than 
Eqn. 1361 due to the oscillatory terms in the integrand 
canceling out behavior in the exponential tail. This solu- 
tion encompasses the long time behavior (t ^> íl _1 ), and 
breaks down if p sym is singular for oj > Q, or if u < C 46 . 
The lack of dependence of the effective Rabi frequency on 
the detuning follows naturally from situations where the 
density of states at —8, p(—6) is sizeable: that portion of 
the density of states is resonant, and the resonant behav- 
ior dominates. When u < 0, we may be in the detuning 
dominated regime, and these results no longer hold. 



C. Discussion 

The long-time behavior, Eqn. 1401 may seem peculiar at 
first. The power-law decay envelope comes from the por- 
tion of the bath density of states that is "on-resonance" 
with the oscillations. In particular, for short times, Rabi 
oscillations are insensitive to detunings less than the Rabi 
frequency. This window narrows for longer time, leading 
to the observed power law behavior. More curiously, for 
times t> m, there is an overall phase shift of ir/4 that is 
independent of many of the parameters of the bath. This 
dynamic shift eludes a similar immediate qualitative de- 
scription, arising from the continuity of the bath density 



of states and the pole induced by probing of the bath (by 
means of the applied Rabi field) that is insensitive to first 
derivatives of the bath, i.e., bath variables come in only 
in their quadrature and higher moments. 

It is crucial to distinguish these results from more Stan- 
dard inhomogeneous broadening results. Unlike the case 
of Doppler broadening, for example, in this system the 
individual system detunings do not change as a func- 
tion of time. When the correlation function is short- 
lived, the behavior would instead follow directly from the 
well undcrstood results describing inhomogenous broad- 
ening duc to Doppler shifts of lascr-induccd Rabi oscil- 
lations. To extend these results to describing inhomoge- 
neous broadening in ensembles of self-assembled quan- 
tum dots, for example, would require considering the ad- 
ditional effect of inhomogeneity in the Rabi strength, Í2, 
and it is unclear that the corresponding phase shift would 
behave in a similar manner to the case of a single Rabi 
frequency. As such, that analysis is beyond the scope of 
this paper. 

To compare these analytical results to exact solutions, 
we solve the system for finite spin, where we exactly eval- 
uate the trace for finite N. In the homogeneous case, the 
Dicke picture of collective angular momentum allows us 
to go to large N, while the inhomogeneous case requires 
exponential operations for the exact value, but can be 
simulated through stochastic modeling of the trace func- 
tion. Wc compare several different vàlues of inhomogene- 
ity and bath size N in Figure""!] and show strong quan- 
titative agreement between the short time and long time 
approximations with the numerical evolution. We also 
show the convergence to full contrast for a fixed time 
and varying Rabi power. 

The previous analysis assumed that the bath density 
of states could be taken as a Gaussian, which can be 
shown to reduce to the assumption of Gaussian noise 4 ^. 
However, inclusion of non-Gaussian effects leads in the 
general cases to integrals without easy solution. 

The deviation from Gaussian noise for the fourth order 
at T 



00 IS 



3(À 2 Z ) 2 



-2/N. 



(43) 



To investigate the role non-Gaussian effects play, sim- 
ulations with finite numbers of spin were conducted; 
the thermodynamic limit recovers Gaussian statistics. 
For small inhomogeneity, Fig. [2t indicates that increas- 
ing spin number decreases the amplitude of fluctuations 
from the expected Rabi signal. Increasing inhomogeneity 
leads to an apparent additional reduction of fluctuations, 
though quantifying this effect remains dimcult. 



D. Inclusion of the Markovian environment 

For completeness, we now add the coupling of the qubit 
to an additional, Markovian environment, with a corre- 
sponding "radiative" decay to the ground state with a 



(a) 




2 4 

t [27T/Í1] 




2 4 

t [27T/Í1] 



(c) 



04 




2 4 6 

t [2tt/í1] 



FIG. 4: (a) Comparison of short time and long time the- 
oretical results to analytic results, with a a — a/2, X = 1, 
N — 20, and = 2ir. The exact solution (sòlid line) com- 
pares well with the short time decay (dashed line) and the 
long time tail plus phase shift (dotted line). (b) Comparison 
of the same exact result with the case with 7 = 0.5 (dot- 
ted line). (c) Fluctuations for N =6 (sòlid), 10 (dashed), 14 
(dashed-dotted), and 20 (dotted). 



rate 7. The optical Bloch equations must now be modi- 
fied to include this decay, and exact solutions, while avail- 
able, are cumbersome. However, the case of FID, with 
íl = 0, is immediately understandable. In the Heisenberg 
picture, we write 



dS+ 
~dT 



-7S+ + i\A z S + + jF + , 



(44) 



where the stochastic input field !F+ has the Standard 
Markovian kernel, (P-(t)P + (t')) = 2jS{t - £')■ Solving 
this equation of motion exactly for (S+(t — 0)) = 1/2 



gives 



FID,~f 



(S+(0)> 



exp(— i\A z t — jt) 



bath 



The Markovian environment results in irreversible expo- 
nential decay, which would be unaffected by any spin- 
echo experiment. For long correlation time baths, spin- 
echo experiments provide a direct measure of the Marko- 
vian component of the decay. 

For the case with Í7 ^ 0, we use the analytical evolu- 
tion for solving the Kraus operator form for the Bloch 

vector, dS/dt = MS + u and evaluating the expectation 
vàlues of all functions of A z at each time point. The 
results are cumbersome, and checked only numerically. 
Qualitatively comparing the result, shown in Figure^J), 
to the 7 = case, we see that the crucial features of phase 
shift and fast initial damping with long-lived oscillations 
remain, but that the steady state population difference 
is shifted, and the oscillations are eventually damped by 
the exponential tail. The corresponding lineshape is ap- 
proximately 
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O7 



[( 7 /2) 2 + (xa z + sy + 02] 



(46) 



Finally, we now consider the effect decorrelation of A z 
plays in the observed oscillations. Unlike the case of free 
evolution, [Ho(t), Ho(t')] 7^ 0, leading to difficulties in 
analytical evaluation of the propagator including decor- 
relation effects. We can, however, consider it partially 
by use of a Magnus expansion, in the limit of weak to 
intermediate strength mesoscopic spin bath. 

Formally, we transform Hq to the interaction picture 
with respect to ClS x , and with 5 = 0, giving 



H = XÀ z [cos(ÍU)S z + sm(nt)Sy 



(47) 



If we assume A z varies on time scales much longer than 
1 / O, we can use a magnus expansion. The expansion 
of this Hamiltonian at peaks in the expected oscillations 
(t = 2tt/í7) ycilds 



with 



U ef f (r) = exp[-zr(ií + Hi + . . .)] (48) 



H = (49) 
H x = 2(XÀ z ) 2 /nS y . (50) 

We neglect higher order contributions by assuming that 
(Al) 4 X 4 /n 4 -C 1. Thus, the envelope of oscillations 
should be determined by 



u eff(t) = exp[-z 



2(XÀ z (t')f 

n 



dt S.,, 



(51) 
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FIG. 5: Expected value of 2S Z as a function of time for 
an initial state S z = 1/2, with driving and low frequency 
noise A = fi/2. Decorrelation rates of (black), fi/100 (red), 
fi/30 (green) and fi/10 (blue) are shown. It is apparent that 
decorrelation leads to steady state vàlues of zero, faster than 
1/vt decay, and reduced phase shift (from 7r/4). 



where we have assumed that over a time r, A z (t) is fixed. 
We note that this form of quadratic noise has been in- 
vcstigatcd in dctail 20 . The lowest order perturbative ex- 
pansion requires determining the spectral function 



Sa 2 M 



dr 
2^ 



(À z (t + r) 2 À z (t) 2 )e-^ T . (52) 



Assuming Gaussian statistics for À z , we find 

(Az(t) 2 À z (t') 2 ) = ( J duS(üi)f + 2(J dujSiíüy^-*")) 2 

(53) 

The resulting decay of oscillations should be given by 

sin 2 (wí/2) 



4A 4 



(54) 



Oscillations with low frequency noise for S(SÏ) are shown 
inFig.0 

We remark that for time independent A z and a Gaus- 
sian density of states we can evaluate the expected enve- 
lope of oscillations, 

(U eff (t)) - (1 + 2itp"(0)/p(0)Ü)- 1 / 2 , (55) 

reproducing the l/yt long-time tail found by more exact 
means. Furthermore, as (A'l) ^ in general, a phase 
shift due to small rotations around S y after each Rabi 
oscillation should be expected. 



V. APPLICATION TO EXPERIMENTAL 
SYSTEMS 

We now consider the physical systems where the quasi- 
static bath assumptions are vàlid. One case is nuclear 
spins interacting with an electron spin in a quantum 
dot. There, the dominant internal bath Hamiltonian is 



aligned with the bath-interaction, and [Hg, A z ] — iden- 
ticall y 1 i. Another system is superconducting qubits, 
where charge-traps and other few-state fiuctuators play a 
similar role. Finally, some NMR systems with one probed 
spin (e.g., a Carbon-13) is coupled to many spins (e.g., 
nearby hydrogens). 



A. Electron spin in quantum dots 

The hyperfine interaction between the electron and the 
nuclear spins is given by 

Hqd = J e B S z + Xqd ■ ï+ "í" B ° J2 % ( 56 ) 



where the oc |?/>(rfc)| 2 are weighted by the 

electron wavefunction's overlap with lattice site k. 
For GaAs quantum dots with N nuclear spins, 
Xqd — 207ns~ 1 /v / ï\^*L, with the normalization condi- 
tion J2k a 1 = 1- Idcntifying 7 e i?o = ui and adding a 
time dependent ESR field in the x-axis the model Hamil- 
tonian is recovered, with A = Xqd- In the case of a large 
applied magnètic held (us ^> Xqd), the results heretofore 
derived hold. 

Whilc experiments investigating ESR and spin dephas- 
ing in single quantum dots are ongoing, the contraints 
imposed by the quasi-static limit make Standard ESR 
technically dimcult to achieve. For small quantum dots 
(e.g. single electron quantum dots), Xqd > 0.1 ns _1 
(N ~ 10 5 ). If the decorrelation rate T is determined by 
dipolar-diffusion processes, it can be no faster than 10 
ms^ 1 , the linewidth of NMR in bulk GaAs, and the quasi- 
static limit is appropriate. For the quasi-static held to be 
overcome, either many measurements must be taken in 
the correlation time of the bath, or ESR held's Rabi fre- 
quency Í2 must be much larger than 0.1 ns . This latter 
requirement may be quite dimcult for current experimen- 
tal parameters. Active correction sequences may lower 
this rate, by averaging through NMR pulses the dipolar 
interactioni&íííi, making it possible to perforin accurate 
phase estimation within the enhanced correlation time of 
the bath. 

There is another important situation where the meso- 
scopic bath picture holds: during an exchange gateü In 
particular, for two tunnel coupled quantum dots, each in 
the single-electron regime, the overall Hamiltonian is 



H 



DQD 



h qd + Hqd + JS 1 ■ S 2 



(57) 



where J is the exchange interaction between the two dots. 
While the m s — ±1 triplet states (| fí), | ü)) are well 
separated by Zeeman energy, the m s = states (| tl), | J. 
T)) are nearly degenerate. In the basis {| tl), | IT)}, this 
two-level system is described by the Hamiltonian 



H. 



-XA 



m s =0 



(DQD) 

n/2 



j2 ü/2 



XA 



(DQD) 



/2 



(58) 



where \À Z DQD) = E fc , i=1 , 2 A^RX^ is the bath 
parameter for this two-level system, and fi = J. This 
Hamiltonian maps exactly to Equation|SJ 

Recently experiments demonstrated coherent oscilla- 
tions driven by controlled exchange interactions 8 . The 
system was prepared in an eigenstate of the mesoscopic 
spin bath interaction (e.g., the | fj.) state) and coherent 
oscillations between | fj.) and | JJ) were driven by applying 
a pulse of non-zero f2 for finite time. These experiments 
correspond exactly to the driven evolution examined in 
the previous section. As such, for small exchange vàlues 
(fi ~ A) wc expect the exchange operations to exhibit 
the characteristic power-law decay and phase shift found 
in this work. Our approach agrees with system-specific 
theoretical predictions of Ref. 129 . 



B. Superconducting qubits 

The quasi-static assumptions also hold for bistable 
two-level fluctuators, where the I z eigenvalues denote 
which state is populated, and the dominant interaction 
with a charge system is a conditional capacitive energy, 
leading to a S Z I Z interaction between the two. Fur- 
thermore, l//-type two-level fluctuators, considered the 
source of most Josephson junction dephasing^iiS, are 
quasi-static in the high frequeney regime, i.e., above the 
bath cutoff. Curiously, the spin-boson modeiS does not 
fulfill this criterion. Other pseudo-spin systems, such as 
flux qubits, may or may not have this structure, depend- 
ing on the bath degrees of freedom. 

The extension of our results to superconducting qubit 
experiments is motivated by results characterized by low 
contrast, long-lived Rabi oscillations, similar to those ex- 
hibited in our model. For a system with some minimum 
time of applied Rabi field (as is often the case in these ex- 
periments) the initial, fast decay could be encompassed 
by oscillations oceurring within this minimum time, while 
the observable oscillations correspond to the long-time 
tail. 

Of the superconducting qubit experiments, we focus 
on Ref. E as 1/f noise is a dominant term for simi- 
lar devices in CW behavior. In addition, there are ob- 
served resonances with flucuators in the measurement 
spectrumià, which suggests that coherent coupling to 
these two-level systems may be possible. A model which 
includes preparation and measurement errors is consid- 
ered in Appendix [Bj we note the results here. We find 
that in the low power regime, our model explains the ob- 
served behavior, but at higher powers incoherent heating 
due to driving microwaves play an important role as well. 

Fitting the parameters from Appendix IBI to the Rabi 
oscillations shown in RefJ£ yield the vàlues in TableQ] 
The results are consistent with bath domination for low 
power (fi < 0.5ns _1 ) and heating / preparation error 
dominating at high power. This gives the characteristic 
low contrast oscillations observed, even though the static 
measurement emeieney is quite high (order 75-85%). The 



TABLE I: Fitted vàlues (explained in Appendix [ÏÏJ for the 
model describing the experiments of Ref. Usi The vàlues in 
parantesis are the 90% confidence intervals of the fit. 

|Mt| 0.75(4) f 

Ml 1.00(3) 

T 0.10(4) ns" 1 

A 0.27(1) ns _1 
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FIG. 6: Plot of oscillations reproduced from Ref. [13 (filled 

circles). The probability of spin flip, P\ — ~ — {S z (t)) is 

1/2 

plotted versus Rabi frequeney íl oc Pjq a fixed Rabi evolution 
time of 25 ns. Overlaid for the given fit parameters are the 
expected vàlues from our model (sòlid line). 

resulting oscillations are compared to Ref. [T3 in Fig. |ü| 



C. Other systems 

Under certain conditions, NMR systems may be de- 
scribed by a quasi-static mesoscopic spin bath . For 
example, large molècules or crystalline structures with 
long correlation times and coherent dynamics within the 
molecule, where one or few spins are coupled to many ad- 
jacent spins in a static configuration, may be described 
by our model. In contrast, liquid state NMR with small 
molècules has bath characteristics dominated by tran- 
sient couplings to other molècules, and our model is in- 
appropriate. 

One example is in cross-polarization experiments on 
solid-state ferrocene, where defocusing sequences are ap- 
plied to generate a Hamiltonian equivalent to Eqn. [F] 54 . 
By direct observation of the cross-polarization signal, 
driven evolution may be observed. We remark that previ- 
ous experiments, investigating spin-bath coupling in frec 
evolution 43 will not show this behavior. System specific 
workSà is in general agreement with our results. 



VI. METHODS TO MITIGATE BATH EFFECTS 

While the strongly coupled mesoscopic spin bath leads 
to dephasing during free evolution and reduced contrast 



and a dynamical phase shift during driven evolution, the 
long-time correlation of the bath allows for correction. 
We start by considering spin-echo-type passive error cor- 
rection, which can be implemented if a sufficiently strong 
Rabi frequency for driving is available. We then considcr 
how entropy transfer from the quasi-static bath to the 
cnvironment (bath-cooling) may be an effective alterna- 
tive means of reducing the damping due to the bath. In 
essence, for specific bath models, it may be possible to 
prepare coherences in the bath which effectively reduce 
the uncertainties engendered by the bath. The inherently 
quantum mcchanical nature of this leads us to considcr 
this mechanism as environmental squeezing. 

A. Passive error correction 

With strong Rabi power (fi ^> A), high fidelity Rabi 
oscillations are possible even in the presence of the meso- 
scopic spin bath. This is the power- broadened limit. The 
error of such pulses scales as (A/fi) 2 and as such deter- 
mining the quasi-static bath strength for a given system 
will yield the required Rabi power to overcome the effects 
of the bath through passive error correction. 

In this limit, spin-echo techniques work well to greatly 
extend the lifetime of oscillations, and passive error cor- 
rection 31 will be straightforward to implcmcnt. Higher 
order pulse sequences than those considered in Section III 
may lead to further improvemen1j£&. Other techniques, 
such as quantum bang-bang 3 ^, are already implícit in 
this analysis in the form of coherent averaging theory^à. 

B. Cooling the bath 

When experimental limitations of Rabi power or other 
experimental difficultics makc minimization of the cou- 
pling strength impossible, the quasi-static bath remains a 
problem uncorrectable by passive error correction. Then, 
active cooling techniques may be a useful alternative. We 
now outline how, for the specific case of sealar interac- 
tions with the bath, qubit manipulation may lead to ef- 
icient entropy transfer from the mesoscopic spin bath to 
the qubit system. By exploiting resonances between the 
bath and system, we may cool the bath degrees of free- 
dom. Furthermore, in the quasi-static limit, coherences 
developed between bath spins can lead to additional im- 
provements through a form of squeezing. 

To illustrate cooling, the explicit case of an electron 
spin in a quantum dot is illuminating. For that system, 
the full spin-coupling (including the terms previously ne- 
glectcd in the rotating wave approximation) is a sealar 
interactioré, 

\S ■ À = \S ■ (^Ça k ï>^J . (59) 
The dominant terms in the internal Hamiltonian of the 



nuclear spins is 

ff A = (60) 

fe 

due to the Larmor precession of the nuclear spins in the 
external magnètic field. When the effective splitting of 
the qubit (electron spin) is comparable to u>k, the elec- 
tron spin comes into resonance with the nuclear spins. 
In essence, by tuning the system near such a "noise- 
resonance" , the bath can be driven in a double-resonance 
manner (analogous to electron-nuclear double resonance, 
ENDOR) , where the on-axis (z — z) coupling is offset by 
the changing frequency of the applied Rabi field. Before, 
the internal spin dynamics were considered to be much 
slower, i.e., ujk *C \u) + A(A Z )|. For double-resonance, 
this condition is no longer satisfied, and the quasi-static 
approximations begin to break down. 

For this resonant regime, the reader is referred to pre- 
vious work on electron spin cooling 23,57 for the details 
of the process. The benefits of cooling are outlined 
here. In the fully cooled state, the quasi-static bath 
acts entirely as an additional detuning on the system, 
as (A z ) 2 = (A 2 ). Near such high polarization, the quasi- 
static bath effects do not supress Rabi and Ramsey-type 
effects, as the bath strength parameter at high temper- 

ature {\\J (À 2 )) is greatly reduced at low temperature. 
Quantitatively, this improvement is given by the ratio of 
Jfid (Eqn. |20J at high and low polarizations, which is 
r) = — P) for a polarization P of the bath (ne- 

glecting inhomogeneous corrections). In terms of exper- 
imental observables, it should result in an increase in 
the observed contrast of Rabi oscillations and improved 
Ramsey signals as if the bath strength parameter, A, had 
been reduce by the improvement factor 77. The improve- 
ments will be limited by coherences developed within the 
bath, by the correlation time of the quasi-static bath, 
which will limit the lifetime of the cooling, and by other, 
non-local (and correspondingly difficult to cool) degrees 
of freedom, implicit in the Markovian cnvironment com- 
ponent of this model. 

Cooling could proceed along the following lines. The 
system is prepared in a | |) state by coupling to the 
Markovian environment, e.g., by moving to a strongly 
coupled region of parameter space. This is the current 
preparation method for superconductor-based qubit de- 
signs; in quantum dots, the energy difference between 
spin states due to zeeman splitting or exchange inter- 
actions allows preparation by controlled coupling to a 
nearby Fermi seaè. Then a weak, negatively detuned 
Rabi field is applied for a short time (order one Rabi 
flop). This process is repeated many times. Including 
state measurement after the Rabi field will allow mea- 
surement of the cooling efficiency, as the bath will lead to 
additional flip to the | f) state in the far detuned regime. 
In some cases, where the natural freqüències of the bath 
degrees of freedom are large, resonances between the bath 
and system may be observed. Tuning to the "red" side 
of these resonances will also force cooling to oceur. 



As cooling is effected only through the logical basis 
states (the qubit states of up and down), it is not neces- 
sarily efficient. For a simple scalar coupling (e.g. a hy- 
pcrfine interaction) , cooling proceeds due to the action 
of the A_ operator, where 

A_=J> fc /* (61) 

k 

However, states \T>) with the property A_|2?) = can- 
not cool under this action. This is similar in practice 
to single mode cooling, and as such it does not greatly 
improve the bath polarization or entropySà for intermedi- 
ate or larger bath sizes. However, these states do have a 
uscful symmetry property, namely they are "dark" under 
the action of the collective lowering operator. Thus, the 
natural state of the density matrix after many cycles of 
cooling is 

p = J2pW\v(P)) (62) 



where the sum is over other degenerades (see^ààà for 
details). We now indicate how cooling limited to po- 
larizations such that rj remains order unity can still be 
sufhcient to greatly reduce the dephasing induced by the 
mesoscopic spin bath. 

In particular, if this cooling can be applied along the 
axis of Rabi oscillations (rather than the z axis), the ef- 
fects of such marginal cooling are immediately apparent. 
In particular, the component of noise in the x and y axes 
are averaged by the bias, u>, while the z axis is only av- 
eraged by the driving field, fi. If the z component of the 
bath noise can be reduced, the corresponding Rabi and 
Ramsey signals will be improved. In the case of homo- 
geneous coupling (pi^ = 1/VÏV) to the environment, it 
can be shown immediately that cooling into dark states 
decreases the quadrature in the non-cooled axes, i.e., 

(Àl) = 1/VN (63) 

which is a y/~N improvement over the uncooled state. As 
the other quadratures correspond to noise suppressed by 
the bias, the effect on experiment is immediate and possi- 
ble quite useful. Such Rabi-axis cooling can be achieved 
for nuclei in quantum dots by first cooling along the nat- 
ural (z) axis, followed by a NMR 7r/2-pulse to rotate the 
nuclear coherences to be parallel to the driving field. 

The coherences between spins in the bath have a lim- 
ited lifctime, determined by the decay of the correlation 
functions of the bath. Therefore, for this process to be 
useful, the qubit-bath interaction must be much faster 
than the bath decorrelation time, but this is exactly the 
quasi-static approximation already made. Before each 
run of a Rabi oscillation-type experiment, a cooling se- 
quence could be applied, to keep the bath degrees of 
freedom. The resulting oscillations should reflect the 
"cooled" bath results, and as such could provide a sub- 
stantial improvement over uncooled systems. 



VII. CONCLUSIONS 

The response of a two-level system (a qubit) with 
driving fields coupled to a local, mesoscopic collection 
of finite-level systems and with additional coupling to 
a larger, Markovian environment has analyzed in this 
paper. The natural simplifications of the Hamiltonian 
lead to a quasi-static regime, wherein bath correlation 
functions are long-lived when compared to experimental 
manipulation and measurement time seales. However, 
decorrelation of the bath and dynamic thermalization 
processes lead to apparent decay when a single-qubit is 
manipulated and measured repeatedly. The large vari- 
ations of bias produced by the quasi-static bath can be 
difficult to correct with limited Rabi frequeney ranges 
accessible to experiments. 

The simple case of free induction decay demonstrates 
the role increasing bath strength and non-Markovian ef- 
fects have to play on the short term, non-exponential 
decay, while finite number effects can lead to revivals 
and fluctuations, distinctly mesoscopic effects. This il- 
lustrates explicitly the non-Gaussian nature of the bath. 

In the presence of weak fields, a line-shape distinctive 
to the spin-bath model is discovered, though inclusion of 
additional (Markovian) decoherence and dephasing re- 
covers the expected Lorentzian behavior. As the driving 
field's strength is increased, Rabi oscillations are possi- 
ble, but for short times there may be a phase shift of 
order 7r/4 and fast initial decay of the oscillation enve- 
lope. However, longer times, with only a power law decay 
going as í -1 / 2 , will most likely be dominated by other, 
exponential decay processes. As a result, a quasi-static 
bath could lead to reduced contrast of oscillations that 
may be consistent with experimental observations of su- 
perconducting qubit devices. In addition, experiments 
with single- and double-quantum dots may probe this 
quasi-static bath model more directly. 

Finally, the role of passive error correction, quantum 
bang-bang, and environment cooling have been consid- 
ered, with each appropriate for different ranges of pa- 
rameters. It may be that as the understanding of quasi- 
static, local baths increases the means to mitigate their 
effect will become more apparent. This is a first and 
crucial step towards that end. 
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APPENDIX A: INTEGRALS INVOLVED IN 
DRIVEN EVOLUTION 

We now solve the integrals necessary to evaluate (/) 
and (C(í)) for a given bath density matrix, pbath- Tak- 
ing the continuum limit for this density matrix, we can 



evaluate a function h(XÀ z ) as 

Tr[h{\À z )p bath ] = £/i(A)^ 



(Al) 



dAh(A)p(A) (A2) 

Accordingly, the long time population difference, given 
by (/) has the form 

1 



(/) = 1-E^<A| 



1 + (S±M*)2 



|A) (A3) 



where A is the eigenvalue of XÀ Z . Going to the continuum 
limit, and assuming the power is weak (appropriate for 
lincshape measurements) , 

r°° n 2 
(/> = 1 -L dAp(A) (AWTí? (A4) 

Furthcrmore, the integral depends only upon Cj — 
y/(A + S) 2 + íl 2 , so we define the symmetric density, 



p sym (u) = p(-6 - V^ 2 - n 2 ) + p{-s + ^Juj 2 - n 2 ). 



Then, 



(/>=1 
l-íl 1 ' 2 



du- 



íl 2 



in 

dd 



vw 



2 -n 2 ' 



(w/n + i)V(w/n + 2)(s)' 



(A5) 

(A6) 
(A7) 



While this integral can be solved for certain density of 
states, we note that for small freqüències (< n) a rel- 
atively fiat density of states can be approximated in a 
Taylor series, while for large freqüències, the 1/Cj 2 behav- 
ior kills the higher freqüències components of the bath. 
Expanding all of these, except for the pole at Cj = 
the non-singular, non-oscillatory terms in Ecin. lAíU near 
Cj = are 



Psym(Cj + n) 



(w/í2 + l) v /w/í2 + 2 



= V2p{-5) 



1 — UUJ - 



where the time scale is 



5 

3n 



: V2p(— S) exp [— uQ] 
Üp"{-8) 



PÍ-S) ' 
Solving the integral yields 

(f)~l-p(-6)y/2TrCl/u 



(A8) 
(A9) 

(A10) 
(AH) 



The oscillations at large t can be evaluated in a manner 
similar to the lincshape. 



(C(í)> 



du — , 
n Vw 2 - n 2 



(A + 5) 2 + n 2 

n 2 exp(— iojt) 



psym(w)- 



(A13) 



n x /2 e 



o w (w/n + i)^/(w/n + 2)(w)^ 



where we have transformed once more to an offset vari- 
able, üj — lü — n. For this final integral we use the station- 
ary phase approximation at long times, requiring that 
Psym(y> > n) not be singular. 

This linear term (Eqn. IA9|) determines the corrections 
to the stationary phase integral. We now define a time 
dependent angle and an effective time, 



9{t) = 1/2 tan" 1 [t/u] 
r(í) = Vt 2 + u 2 . 



(A15) 
(A16) 



Using these substitutions, we can solve for the stationary 
phase result exactly: 



(C(*)> 



2ílp{-6) e - int 
= p(-S)e- mt 



duj- 



(A17) 
(A18) 



APPENDIX B: COMPARISON WITH 
SUPERCONDUCTING QUBITS: MODEL 

We develop a theory describing supcrconducting qubit 
systems that includes both the spin-bath model of this 
paper and incoherent heating from microwave pulses used 
to implement driving. To characterize the experiment, 
a simple model involving several static parameters was 
used. One is the preparation efnciency, /, which is the 
probability of preparing "down" (the desired initial state) 
minus the probability of preparing "up" (through bad 
preparation, probably due to high temperature). I is 
set entirely by the microwave temperature. Assuming 
that the only signficant source of heating is the incident 
microwave field, Stefan's law gives T oc Accordingly, 
/ = tanh(yïyn), where T 2 is a measure of the cooling 
power. 

In addition, there is a finite probability of measuring 
the wrong state, given by the conditional probabilitics, 
-Mfj and M^. The first is the probability of measur- 
ing up if the system is in the up state at the time of 
measurement, and the second is the same, but for down. 
Assuming that Jl4fj + JVÍfj_ = 1 and similar ly for tíiej. case, 
the measured signal is then 

2s z = A% - M u + / p S4)(2A% + 2M Ü - 2) (Bl) 

where S z (l) is the result given by the model of this paper 
for an initial j state (sec Eans. lATTl and IÀ18I and recali 

that 25 z = Re[(/ + C(í))])- 

To include bath effects, a Gaussian density of states 
(appropriate if the size of the bath is > 20 effective spins) 
is assumed. Then p(8) = exp[— <5 2 /(2A 2 )]/\/2^A 2 , where 
A is the bath strength parameter. While non-Gaussian 
effects, leading to mesoscopic fluctuations, are important, 
their inclusion leads to too many highly correlated fitting 
parameters. 
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